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ABSTRACT 

Strong tidal interaction with the central star can circularize the orbits of close-in planets. With 
the standard tidal quality factor Q of our solar system, estimated circularization times for close-in 
extrasolar planets are typically shorter than the ages of the host stars. While most extrasolar planets 
with orbital radii a < 0.1 AU indeed have circular orbits, some close-in planets with substantial 
orbital eccentricities have recently been discovered. This new class of eccentric close-in planets implies 
that either their tidal Q factor is considerably higher, or circularization is prevented by an external 
perturbation. Here we constrain the tidal Q factor for transiting extrasolar planets by comparing their 
circularization times with accurately determined stellar ages. Using estimated secular perturbation 
timescales, we also provide constraints on the properties of hypothetical second planets exterior to 
the known ones. 

Subject headings: planetary systems 
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1. INTRODUCTION 

The median eccentricity of the current sample of ~ 300 
planets is 0.19, while it is 0.013 for close-in planets with 
semi-major axis a < 0.1 AU. The circular orbits of close- 
in planets most li kely result from o r bital circularizatio n 
due to tides (e.g., iRasio et alj|1996t iMarcv et al.lll997l ). 
This requires the tidal circularization time r c i rc to be 
short compared to the age of the system r agc . Since 
r C i rc is a very steep function of a (see Eq. Q] or [6| , while 
r agc ~ 1 — 10 Gyr for most systems, a sharp decline in 
eccentricity is expected below some critical value of a. 
However, the observed transition seems to occur around 
0.03 — 0.04 AU, whereas the calculated T c i rc becomes com- 
parable to r ago at ~ 0.1 AU as we see below (Figured]). 
Since almost one quarter (presently 16/68) of planets 
within 0.1 AU have e > 0.1, their high eccentricities de- 
mand explanation. 

First, we calculate the circularization times for transit- 
ing planets, and compare them with the estimated ages 
of the systems. The circularization time r C i rc = — e/e, 
where e is the sum of the eccentricity change due to the 
tides raised on the star by the planet and those raised on 
the planet by the star, is dGoldreich fc Soterlll966t iHutl 
llMfl lEggleton et al.lfl998t IMardling fc LiiJl2002j ): 
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The subscripts p and * represent the planet and star, re- 
spectively. The modified tidal quality factor for a planet 



is defined as Q 1 



3<5p/2fc p , where fe p is the Love num- 



ber, and Q p is the specific dissipation function, which 
depends on the planetary structure as well as the fre- 
quency and amplitude of tides. We also define 
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Generally F* and F p are comparable, and thus the stel- 
lar damping is negligible unless the planet-to-star mass 
(radius) ratio is large (small) or <C Q' p . We define the 
circularization time due to damping in the planet as 
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where n = y/G(M* + M p )/a 3 is the mean motion, fi rot 



Note that T c j rc .o can be shorter or longer than r c ; rc , 
depending on the sign of which changes at 

(fi*,rot/rc)crit = 18/ll(/i// 2 ). In the limit e -> 0, this 
equation leads to the stand ard expression for the circ u- 
larization time (Eq. 4.198 of lMurrav fc DermottH"l999D . 

Figure [T] compares the circularization times calculated 
from Eq. [T] and [S] with the estimated stellar ages for 
the systems in Table [T] 1 . Here we assume Q' p — 

10 5 , and Q'x — 10 6 , which are the standard values 
motivated by measurements in our Solar Sys tem (e.g. 
lYoder fc PeaIi[T98l iZhang fc HamiltoiJ 12008). an d for 
main-sequence stars (e.g. ICarone fc Patz old 2007). For 
^p,rot, we assume that planets with circular orbits are 
perfectly synchronized, i.e., fl p ^ot/n = 1, si nce the spin- 
orbit synchronization times are ~ 10 -3 T c j rc (Rasi o et al.l 
1996). On the other hand, planets with eccentric orbits 
should spin down until th ey reach quasi-synchronization 
(|Dobbs-Dixon et aTll2004[) : in practice we adopt a plan- 
etary spin frequency such th at the rate of change o f spin 
frequency is zero (Eq. 54 of IMardling "fc~ Lin 2002]) . For 

1 Note that, in Eq. [T] it is implicitly assumed that the star 
and the planet both have zero obliquity. Currently available mea- 
surements of the Rossiter— MacLaughlin effect show that the plan- 
etary orbits in general are closely alig ned with the stellar equator 
HQueloz et al. 2000: Winn et al. 2005). Current exceptions may be 
the H P 17156 and XO-3 systems iNarita et all2007l;rHebr'ard et al.l 
2008), which we exclude from our analysis. 
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Fig. 1. — Circularization times calculated from Eq. \E\ (stars) , 
Eq. [l]with f2* >ro t/n = 3d (orange circles) and 70 d (blue circles), 
compared with the estimated stellar ages (squares) for systems with 
transiting planets. Here we assume Q' p = 10 5 and Q' t = 10 6 . For 
all transiting planets, the estimated circularization time is shorter 
than the age of the host star. 

the stellar spin, we assume typical peri ods derived f rom 
the observed v TOt sini, P* )ro t ~ 3-70d (|Barnes!l2001l ). 

Figure [T] compares the circularization times calculated 
from Eq. [1] and [HI with the estimated stellar ages for the 
systems in Table [TJ For most systems the tidal damping 
in the star is negligible. Two systems with non-negligible 
stellar damping are WASP-14 and HAT-P-2. Both have 
a large planetary mass (see Table [1} , and thus the first 
term in Eq. [1] is significant. Clearly, the estimated cir- 
cularization times are always shorter than stellar ages, 
which implies that all these planets should have been 
circularized by now if their tidal Q values were similar 
to those of their solar-system analogues. However, our 
sample contains at least 6 systems with a non-zero or- 
bital eccentricity. Possible explanations are that (1) the 
structure of these planets is different and their actual Q 
is greater than what we assumed; or (2) they currently 
experience external perturbations which maintain their 
orbital eccentricities against tidal dissipation. 

2. CONSTRAINTS ON THE TIDAL Q VALUE OF CLOSE-IN 
PLANETS 

Now we use Eq. Q] to place constraints on the tidal Q 
values of the planets. An upper limit on Q is provided 
for planets with circular orbits since the circularization 
must have occurred within the lifetime of the systems 
( r circ < T ago)- Our assumption here is that these close-in 
planets formed th rough tidal circu l arizat ion of initially 
eccentric orbits. iNagasawa et al.l (|2008|) showed that 
about one third of multiple planetary systems could form 
close-in planets through tidal circularization following a 
large eccentricity gain through planet-planet scattering 
or Kozai-type perturbations. Direct observational evi- 
dence for initially large orbital eccentricities comes from 
the absence of planeta ry orbits within twice the Roche 
limit around the star (jFaber et al.l 120051 iFord fe Rasiol 
[2001 . 

On the other hand, close-in, eccentric planets impose 
a lower limit on Q values, since r c i rc > r age is expected 
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Fig. 2. — Estimated tidal Q factors for the case of slowly ro- 
tating stars (f2*/n < (f2*/n) cr i t ). Upper/lower limits calculated 
from Eq.[6]are shown in black down/up triangles for planets with 
zero/non-zero eccentricities. Open triangles are the corresponding 
estimates from Eq. [T] which approach black ones as we take 2, 5, 
and 10 times the minimum stellar tidal Qs (blue, orange, and green 
triangles, respectively) obtained from Eq. [9] 

for these systems, provided that they are not currently 
subject to any eccentricity excitation mechanism. 

For -P*. ro t ~ 3 — 70 d, all planets in Table [1] take 
Q*, m t/n < (£l*,rot/™)crit, and hence F* > and T circ , > 
T c i rc . For planets with zero (non-zero) eccentricity, we 



require t c 



^ 7"circ,0 ^ 7"a 



(r. 



■,o)- In 



other words, we assume that zero (non-zero) eccentric- 
ity planets have been (have not been) circularized within 
the lifetime of the system, independent of the rotation 
period of the star. This gives the upper and lower limits 
for circular and eccentric planets, respectively, as 
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The latter also gives the lower limit for the stellar tidal 
Q factor, since the denominator must be positive: 



Q* > f n 



Mp 



R, 



(9) 



This corresponds to a minimum stellar Q value of Q' t ~ 
3 x 10 4 - 4 x 10 7 , with a median value of 0.4 - 1 x 10 6 
for P*,rot = 3 — 70 d, which a gree well with observations 
(e.g. iCarone fc Pateoll[200l . 

Figure [2] shows the upper and lower limits for the plan- 
ets' Q values as a function of semi-major axis for circular 
and eccentric orbits, respectively. Since the lower limits 
on Q values for eccentric planets depend on we take 
three different cases of = 2, 5, and 10 <2'„ min as exam- 
ples. Note that, with this definition of Q*, Q' p becomes 
independent of the stellar spin rate. Since circularization 
times are shorter for planets with smaller orbital radii, 
we tend to overestimate the maximum Q values at the 
shortest-period end. All transiting planets appear within 
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the range 10 5 < Q' p < 10 9 . The figure also shows that the 
high eccentricities of some planets (marked with upper 
triangles) can be explained by assuming relatively large 
(Q'p £ 10 6 ) but reasonable (Q' p < 10 9 ) tidal Q values. 

Although these estimated Q values are larger than 
those of Jupiter or Neptune, they cannot be excluded. 
Recent theoretical studies of the excitation and dissipa- 
tion of dynamical tides within rotating giant planets have 
shown that tidal Q values fluctuate strongly depending 
on the tidal forcing frequency, and the effective Q's could 
go up to ~ 10 9 depending on the spin rate and internal 
structure of the planet (e.g., presence /abs ence of a core, 
radia t ive envelo pe, or a density jump, see lQgilvie fc Lir] 
12004 iWul 120051 ). According to these recent models, it 
appears possible that some planets maintain large eccen- 
tricities simply because of their larger Q values. 

3. DYNAMICAL PERTURBATIONS 

Candidate perturbation mechanisms that could excite 
and maintain planetary eccentriciti es include (1) tidal 
inter action with the central star (Do bbs-Dixon etldl 
2004), (2) quadrupole or higher-order secular perturba- 
tion from an additional body, or (3) resonant interaction 
with another planet. Here we discuss the effects of these 
competing mechanisms against tidal eccentricity damp- 
ing. 

Tidal dissipation inside the central star can increase 
the planet eccentricity only when f2* jro t > n, or equiv- 
alently de/dt > in Eq. [TJ For a synchronized planet 
(^p rot ~ n) with a small eccentricity (e 2 <C 1), we obtain 
Q*, r ' ot /n > 18/11(1 + 7 /l8(Q'JQ' p )(MjM p ) 2 (R p /R*) 5 ). 
For a Jupiter-like planet around a main-sequence star, 
this yields f2* !ro t > 9.6n for ~ Q' p . Since the rota- 
tion period for planet-hosting stars typically lies in the 
range 3 - 70 d, eccentricity excitation may occur for plan- 
ets only if their orbital periods are within 29 - 673 d or 
longer. For a 10Mj planet with radius lRj, we have 
^*,mt > L7n, which corresponds to an orbital period 
greater than 5.1 d. Therefore, this is unlikely to be re- 
sponsible for the eccentricity of observed planets within 
a ~ 0.06-0.18 AU. 

Another possibility is an undetected additional planet 
exciting the eccentricity of the detected planet. If there 
is a large mutual inclination angle (i > 40°) between 
the two pl anets, Kozai -type perturbations can become 
important (|Kozail ll9621. Such highly non-coplanar or- 
bits could result from planet-p lanet scattering a f ter dis - 
sipation of the gaseou s disk. IChatteriee et al.1 (|2007f ); 
iNagasawa et al.l (|2008l ) have performed extensive numer- 
ical scattering experiments and showed that the final 
inclination of planets could be as high as 70°, with 
a median of 10-20°. If i < 40°, octupole pertur- 
bations may still moderately excite the eccentricity of 
the close-in planet. The secular interaction timescale 
of a pair of planets with small mutual inclination can 
be derived from the classical Laplace-Lagrange theo ry 
(|Brouwer fc ClemenceHl961l : iMurrav fc DermottJll999f ). 

For this secular perturbation from an additional planet 
to be causing the large eccentricity of the close-in planet, 
it must occur fast enough compared to other perturba- 
tions causing orbital precession. In particular, GR pre- 
cession and tides are important effects that would com- 
pete against the perturbation from the additional body 
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Fig. 3. — Various secular eccentricity excitation timescales for 
the planet G J 436 b caused by a hypothetical planetary (or stel- 
lar) companion G J 436 c with mass M c and semi- major axis a c . 
The solid black lines show the predicted radial velocity amplitudes 
caused by the undetected companion. The dot-dashed line shows 
the threshold right of which the secular interaction between planets 
with i < 40 deg is suppressed by GR precession, while the three 
red dotted lines are the similar thresholds for the Kozai mechanism 
with the assumed orbital eccentricity of the companion e c = 0.01, 
0.5, and 0.9 f r om le f t to right. We use r DD; TKoza i , and tqr as in 
ITakeda et all l l200Sh : IFabrvckv fc Trem ainc (2007). Thresholds for 
T c i rc and T age are also shown for comparison. 

2 . For detail e d dis cussions see iHolman et al.l (|1997| ): 
iKiseleva et all (|1998[ ). 

Figure [3] illustrates the constraints on the mass and 
orbital radius of the hypothetical outer planet in the 
GJ 436 system. These are set by comparing the GR pre- 
cession and secular timescales. Similar results are seen 
for other systems with eccentric close-in planets. Gen- 
erally GR precession occurs faster than any other per- 
turbation mechanism. In order to induce Kozai cycles in 
the inner planet (left of the dotted lines) while not caus- 
ing radial-velocity amplitudes above the detection limit 
of ~ 5ms _1 (below black lines), the mass upper limit 
of the hypothetical planet is ~ lMNeptune- For near- 
coplanar systems even tighter constraints are placed on 
the properties of the secondary planet (left of the dot- 
dashed line). 

However, one caveat is that our diagram only rules 
out the possibility of hypothetical bodies being cur- 
rently responsible for the high eccentricity of GJ436b. 
Also, while Kozai-type perturbations are almost al- 
ways suppressed by GR precession, eccentricity exci- 
tation through secular octupole pertu rbations may be 
occasionally enhanced by GR effects ()Ford et al.l 120001 : 
I Adams fc Laughlinl 1200(f) . For the case of the GJ436 
system, we have numerically tested the effect of a hypo- 
thetical secondary planet c on the eccentricity evolution 
of the inner planet b. We have found that, within the 
detectable radial-velocity limit < 5m/s, planet c can- 
not excite the eccentricity of planet b from 0.01 to the 
observed 0.15, even if its eccentricity is as high as 0.5. 

2 Although stella r and planeta ry rotational distortions cause ad- 
ditional precession HSterndll93Sh . GR precession dominates unless 
the stellar rotation is on the high end. 
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TABLE 1 

Data are from [http : / /exopianet feu/] Ages are computed using the stellar evolution database in ITakeda et al.1 i|2007f) . 

UNLESS MARKED WITH *. MEDIAN VALUES OF THE DERIVED POSTERIOR AGE PROBABILITY DISTRIBUTION FUNCTIONS ARE PRESENTED HERE, 

TOGETHER WITH THE 95% CREDIBLE INTERVALS IN PARENTHESIS. 



Planet ID 


M p [Mj] 


-Rp [Rj] 


a [AU] 


c 


Ah [M Q ] 


R* [Re] 


Age [Gyr] 


OGLE-TR-56 b 


1.29 


1.3 


0.0225 





1.17 


1.32 


92 (0 20-3 00) 


OOTE-TR-1 1 3 h 


1.32 


1.09 


0.0229 


o 


0.78 


0.77 


1 3 28 f1 1 00 - 1 3 921 




0.072 


0.38 


02872 


0.15 


0.452 


0.464 


6 00 fl 00 - 1 00 "l 


OOTE-TR-1 32 h 


1.14 


1.18 


0306 


o 


1.26 


1.34 


96 10 1 2 - 3 841 


HD 189733 h 


1.15 


1.156 


0.0312 


o 


0.8 


0.753 


8 96 f 1 04 - 1 3 72*1 


TrES-2 * 


1.98 


1.22 


0.0367 





0.98 


1 


5.10 (2.40-7.80) 


WASP- 14 b * 


7.725 


1.259 


0.037 


0.095 


1.319 


1.297 


0.75 (0.5-1) 


WASP- 10 b * 


3.06 


1.29 


0.0371 


0.057 


0.71 


0.783 


0.8 (0.6-1) 


HAT-P-3 b * 


0.599 


0.89 


0.03894 





0.936 


0.824 


0.40 (0.10-6.90) 


TrES-l 


0.61 


1.081 


0.0393 





0.87 


0.82 


11.40 (3.20-13.84) 


HAT-P-5 b * 


1.06 


1.26 


0.04075 





1.16 


1.167 


2.60 (0.80-4.40) 


OGLE-TR-10 b 


0.63 


1.26 


0.04162 





1.18 


1.16 


1.20 (0.16-4.64) 


HD 149026 b * 


0.36 


0.71 


0.0432 





1.3 


1.45 


2.00 (1.20-2.80) 


HAT-P-4 b * 


0.68 


1.27 


0.0446 





1.26 


1.59 


4.20 (3.60-6.80) 


HD 209458 b 


0.69 


1.32 


0.045 


0.07 


1.01 


1.12 


2.40 (0.48-4.60) 


OGLE-TR-111 b 


0.53 


1.067 


0.047 





0.82 


0.831 


5.17 (0.17-13.41) 


HAT-P-6 b * 


1.057 


1.33 


0.05235 





1.29 


1.46 


2.30 (1.60-2.80) 


HAT-P-1 b * 


0.524 


1.36 


0.0553 


0.067 


1.133 


1.115 


3.60 


HAT-P-2 b * 


8.64 


0.952 


0.0677 


0.517 


1.298 


1.412 


2.70 (1.30-4.10) 



This result holds for other systems since they have even 
heavier planets. Therefore, we can safely exclude the 
possibility that these planets obtain their current high 
eccentricities through secular perturbation from an un- 
detected outer planet with < MNeptunc if their orbits are 
initially near-circular. 

Yet another possibility is a res onant perturbation from 
an undetected planet. Recently iRibas et al.l (|2008| ) sug- 
gested that the eccentricity of GJ 436 b might be caused 
by a mean-motion resonance (MMR) with an unseen 
super-Ea rth, but there is little o bservational support for 
this (e.g.H ean fc Seifahrtll2008l ). Also, the combined ef- 
fects of GR precession and MMR are not fully understood 
yet. In any case, it is unlikely that such resonances are 
responsible for all the close-in eccentric planets, consid- 
ering the small fraction of extrasolar multiple planets in 
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4. SUMMARY 

In this letter, we have investigated the origins of close- 
in planets on an eccentric orbit. We place constraints on 
the tidal Q factor of transiting planets by comparing the 
stellar age with the tidal circularization time, and find 
that 10 5 < Q' p < 10 9 , which agrees well with current 
theoretical estimates, can explain these eccentric plan- 
ets. We also show that it is difficult to explain the high 
eccentricities of these planets by invoking a current inter- 
action with an unseen second planet. Our results suggest 
that at least some of the close-in eccentric planets may 
be simply in the process of getting circularized. 
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